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We present the first nucleon-nucleon potential at next-to- 
next-to-next-to-leading order (fourth order) of chiral pertur- 
bation theory. Charge-dependence is included up to next- 
to-leading order of the isospin-violation scheme. The accu- 
racy for the reproduction of the A'^A'^ data below 290 MeV 
lab. energy is comparable to the one of phenomenological 
high-precision potentials. Since A''A'^ potentials of order three 
and less are known to be deficient in quantitative terms, the 
present work shows that the fourth order is necessary and suf- 
ficient for a reliable A'^A'^ potential derived from chiral effective 
Lagrangians. The new potential provides a promising starting 
point for exact few-body calculations and microscopic nuclear 
structure theory (including chiral many-body forces derived 
on the same footing). 

PACS numbers: 21.30-x, 13.75.Cs, 12.39.Fe 



The theory of nuclear forces has a long history. 
Based upon the Yukawa idea [1], first field-theoretic at- 
tempts [2,3] to derive the nucleon-nucleon (NN) inter- 
action focused on pion-exchange, resulting in the NN 
potentials by Gartenhaus [4] and by Signell and Mar- 
shak [5]. However, even qualitatively, these potentials 
barely agreed with empirical information on the nuclear 
force. So, these "pion theories" of the 1950s are gener- 
ally judged as failures — for reasons we understand today: 
pion dynamics is constrained by chiral symmetry, a cru- 
cial point that was unknown in the 1950s. 

Historically, the experimental discovery of heavy 
mesons [6] in the early 1960s saved the situation. The 
one-boson-exchange (OBE) model [7,8] emerged which is 
still the most economical and quantitative phenomenol- 
ogy for describing the nuclear force [9,10]. The weak 
point of this model, however, is the scalar-isoscalar 
"sigma" or "epsilon" boson, for which the empirical ev- 
idence remains controversial. Since this boson is associ- 
ated with the correlated (or resonant) exchange of two 
pions, a vast theoretical effort that occupied more than a 
decade was launched to derive the 27r-exchange contribu- 
tion of the nuclear force, which creates the intermediate 
range attraction. For this, dispersion theory as well as 
field theory were invoked producing the Paris [11,12] and 
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the Bonn [13,8] potentials. 

The nuclear force problem appeared to be solved; how- 
ever, with the discovery of quantum chromo- dynamics 
(QCD), all "meson theories" had to be relegated to mod- 
els and the attempts to derive the nuclear force started 
all over again. 

The problem with a derivation from QCD is that this 
theory is non-perturbative in the low-energy regime char- 
acteristic of nuclear physics, which makes direct solutions 
impossible. Therefore, during the first round of new at- 
tempts, QCD-inspired quark models [14] became popu- 
lar. These models were able to reproduce qualitatively 
some of the gross features of the nuclear force. How- 
ever, on a critical note, it has been pointed out that 
these quark-based approaches were nothing but another 
set of models and, thus, did not represent any fundamen- 
tal progress. Equally well, one may then stay with the 
simpler and much more quantitative meson models. 

A major breakthrough occurred when the concept of 
an effective field theory (EFT) was introduced and ap- 
plied to low-energy QCD. As outlined by Weinberg in a 
seminal paper [15], one has to write down the most gen- 
eral Lagrangian consistent with the assumed symmetry 
principles, particularly the (broken) chiral symmetry of 
QCD. At low energy, the effective degrees of freedom are 
pions and nucleons rather than quarks and gluons; heavy 
mesons and nucleon resonances are "integrated out" . So, 
in a certain sense we are back to the 1950s, except that 
we are smarter by 40 years of experience: broken chi- 
ral symmetry is a crucial constraint that generates and 
controls the dynamics and establishes a clear connection 
with the underlying theory, QCD. 

The chiral effective Lagrangian is given by an infi- 
nite series of terms with increasing number of derivatives 
and/or nucleon fields, with the dependence of each term 
on the pion field prescribed by the rules of broken chiral 
symmetry [16]. Applying this Lagrangian to NN scat- 
tering generates an unlimited number of Feynman dia- 
grams, which may suggest again an untractable problem. 
However, Weinberg showed [16] that a systematic expan- 
sion of the nuclear amplitude exists in terms of (Q/A^)", 
where Q denotes a momentum or pion mass, A^ « 1 GeV 
is the chiral symmetry breaking scale, and > 0. For a 
given order z/, the number of contributing terms is finite 
and calculable; these terms are uniquely defined and the 
prediction at each order is model-independent. By going 
to higher orders, the amplitude can be calculated to any 
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TABLE I. Low-energy constants applied in the N^LO NN 
potential (column ^NN'). The d belong to the dimension- two 
nN Lagrangian and arc in units of GoV^^, while the di are 
associated with the dimension-three Lagrangian and are in 
units of GeV~^. The column '-kN' shows values determined 
from ttN data. 



TABLE II. x^/datum for the reproduction of the 1999 np 
database [38] below 290 MeV by various np potentials. 



NN 



nN 



Cl 
C2 
C3 

_ C4 

di + 
ds 

dz> 
di4 - 



di5 



-0.81 
2.80 

-3.20 
5.40 
3.06 

-3.27 
0.45 

-5.65 



-0.81 ±0.15" 
3.28 ±0.23'' 

-4.69 ± 1.34" 
3.40 ± 0.04" 
3.06 ±0.21'' 

-3.27 ± 0.73'' 
0.45 ± 0.42'' 

-5.65 ±0.41'' 



"Table 1, Fit 1 of Ref . [33]. 



''Table 2, Fit 1 of Ref. [34]. 



desired accuracy. The scheme just outhned has become 
known as chiral perturbation theory (xPT). 

Following the first initiative by Weinberg [16], pio- 
neering work was performed by Ordofiez, Ray, and van 
Kolck [17,18] who constructed a NN potenial in coor- 
dinate space based upon xPT at next-to-next-to-lcading 
order (NNLO; v = 2>). The results were encouraging and 
many researchers [19] became attracted to the now field. 
Kaiser, Brockmann, and Weise [20] presented the first 
model-independent prediction for the NN amplitudes of 
peripheral partial waves at NNLO. Epelbaum et al. [21] 
developed the first momentum-space NN potential at 
NNLO. 

In the 1990s, unrelated, parallel research showed that, 
for conclusive few-body calculations and meaningful mi- 
croscopic nuclear structure predictions, the input NN 
potential must be of the highest precision; i. e., it must 
reproduce the NN data below about 300 MeV lab. en- 
ergy with a /datum « 1. The family of high-precision 
NN potentials [9,22,23,10] was developed which fulfills 
this requirement. Due to the outstanding accuracy of 
these NN potentials, it was possible to pin down cases of 
few-body scattering and of nuclear structure that clearly 
require three-nucleon forces (3NF) for their miscroscopic 
explanation. Famous examples are the Ay puzzle of N-d 
scattering [24] and the ground state of ^"B [25]. 

One important advantage of xPT is that it makes spe- 
cific predictions for many-body forces. For a given or- 
der of xPT, both 2N and 3N forces arc generated on 
the same footing. At next-to- leading order (NLO), all 
3NF cancel [16,26]; however, at NNLO and higher or- 
ders, well-defined, nonvanishing 3NF terms occur. As 
discussed, since 3NF effects are in general very subtle, it 
is only possible to demonstrate their necessity and rele- 
vance when the 2NF is of high precision. 

A''A^ potentials based upon xPT at NNLO are poor in 
quantitative terms; they reproduce the NN data below 
290 MeV lab. energy with a x^/datum of more than 20 



Bin (MeV) 


# of data 


N^LO" 


NNLO*" 


NLO*" 


AV18'= 


0-100 


1058 


1.06 


1.71 


5.20 


0.95 


100-190 


501 


1.08 


12.9 


49.3 


1.10 


190-290 


843 


1.15 


19.2 


68.3 


1.11 


0-290 


2402 


1.10 


10.1 


36.2 


1.04 


"This work. 


''Ref. [35]. 
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which is totally unacceptable. Clearly, there is a strong 
need for more precision, implying that going to higher 
order is necessary. 

It is the purpose of this note to present the first NN 
potential that is based consistently on xPT next-to- 
next-to-ncxt-to-leading order (N^LO; fourth order). We 
will show that, at this order, the accuracy is comparable 
to the one of the high-precision phenomenological poten- 
tials. Thus, the NN potential at N^LO is the first to 
meet the requirements for a reliable input-potential for 
exact few-body and microscopic nuclear structure calcu- 
lations (including chiral 3NF consistent with the chiral 
2NF). 

In xPT, the NN amplitude is uniquely determined 

by two classes of contributions: contact terms and pion- 
exchange diagrams. At N"^LO, there are two contacts 
of order [0(Q°)], seven of ©(Q^), and 15 of O(g^), 
resulting in a total of 24 contact terms, which generate 
24 parameters that are crucial for the fit of the partial 
waves with orbital angular momentum L <2 [27] . 

Now, turning to the pion contributions: At leading or- 
der [LO, 0{Q^), v = 0], there is only the wellknown static 
one-pion exchange (OPE). Two-pion exchange (TPE) 
starts at next-to- leading order (NLO, z/ = 2), and there 
arc further TPE contributions in any higher order. While 
TPE at NNLO was known for a while [17,20,21], TPE at 
N'^LO has been calculated only recently by Kaiser [28]. 
All 2tt exchange contributions up to N'^LO are summa- 
rized in a pedagogical and systematic fashion in Ref. [29] 
where the model-independent results for NN scattering 
in peripheral partial waves are also shown. We use the 
analytic expressions published in Ref. [29]. Finally, there 
is also threc-pion exchange, which shows up for the first 
time at N'^LO (two loops). In Ref. [30], it was demon- 

TABLE 111. x^ /datum for the reproduction of the 1999 pp 

database [38] below 290 MeV by various pp potentials. 



Bin (MeV) 


# of data 


N^LO" 


NNLO'' 


NLO'' 


AV18'= 


0-100 


795 


1.05 


6.66 


57.8 


0.96 


100-190 


411 


1.50 


28.3 


62.0 


1.31 


190-290 


851 


1.93 


66.8 


111.6 


1.82 


0-290 


2057 


1.50 


35.4 


80.1 


1.38 



"This work. 



''See footnote [39]. 
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FIG. 1. np phase parameters below 300 MeV lab. energy 
for partial waves with J < 2. The solid line is the result at 
N^LO. The dotted and dashed lines are the phase shifts at 
NLO and NNLO, respectively, as obtained by Epelbaum et 
al. [35]. The solid dots show the Nijmegen multi-energy np 
phase shift analysis [36], and the open circles are the VPI 
single-energy np analysis SM99 [37]. 



strated that the Stt contributions at this order are negli- 
gible, which is why wc leave them out. 

For an accurate fit of the low-energy pp and np data, 
charge-dependence is important. We include charge- 
dependence up to next-to-leading order of the isospin- 
violation scheme (NL0, in the notation of Ref. [31]). 
Thus, we include the pion mass difference in OPE and the 
Coulomb potential in pp scattering, which takes care of 
the L0 contributions. At order NL0 we have pion mass 
difference in the NLO part of TPE, 717 exchange [32], 
and two charge-dependent contact interactions of order 
which make possible an accurate fit of the three dif- 
ferent scattering lengths, Opp, a„„, and a„p. 

Chiral perturbation theory is a low-momentum expan- 
sion. It is valid only for momenta Q <C « 1 GeV. To 
enforce this, we multiply all expressions (contacts and 
irreducible pion exchanges) with a regulator function. 



exp 



2n 



(1) 



where p and p' denote, respectively, the magnitudes of the 
initial and final nucleon momenta in the center-of-mass 
frame. We use A = 0.5 GeV throughout. The exponent 



FIG. 1, continued. 

2n is chosen to be sufficiently large so that the regulator 
generates powers which are beyond the order {v = 4) at 
which our calculation is conducted; i. e., terms up to 
are not affected. 

The contact terms plus irreducible pion-exchange ex- 
pressions at N^LO, multiplied by the above regulator, 
define the NN potential at N^LO. This potential is ap- 
plied in a Lippmann-Schwinger equation to obtain the 
T-matrix from which phase shifts and NN observables 
arc; calculated. The corresponding homogenous equation 
determines the properties of the two-nucleon bound state 
(deuteron). 

The peripheral partial waves of NN scattering with 
L > 3 are exclusively determined by OPE and TPE be- 
cause the N^LO contacts contribute to L < 2 only. OPE 
and TPE at N^LO depend on the axial-vector coupling 
constant, qa (we use qa = 1-29), the pion decay constant, 
/tt = 92.4 MeV, and eight low-energy constants (LEG) 
that appear in the dimension-two and dimension-three 
ttN Lagrangians (cf. Ref. [29]). In the optimization pro- 
cess, we varied three of them, namely, C2, C3, and C4. We 
found that the other LEG are not very effective in the 
NN system and, therefore, we kept them at the values 
determined from -kN (cf. Table I). The most influential 
constant is C3, which has to be chosen on the low side 
(slightly more than one standard deviation below its TriV 
determination) for an optimal fit of the NN data. Our 
choice for C4, which is substantially above the value de- 
termined in -kN, is necessary to bring the phase shift 
down. 
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The most important set of fit parameters are the ones 
associated with the 24 contact terms that rule the par- 
tial waves with L <2. In addition, we have two charge- 
dependent contacts, which brings the number of contact 
parameters to 26. Since we treated three LEG as semi- 
free, the total number of parameters of the N^LO poten- 
tial is 29. 

In the optimization procedure, we fit first phase shifts, 
and then we refine the fit by minimizing the obtained 
from a direct comparison with the data. The phase shifts 
at N"^LO for np scattering below 300 MeV lab. energy 
are displayed in Fig. 1. The /datum for the fit of the 
np data below 290 MeV is shown in Table II, and the 
corresponding one for pp is given in Table III. The x^ 
tables demonstrate a dramatic improvement of the NN 
interaction order by order. It is clearly revealed that, at 
NLO and NNLO, the reproduction of the NN data is of 
unacceptably poor quality. However, at N^LO, the quan- 
titative character is comparable to the phenomenological 
high-precision Argonne Vis potential [22]. 

In conclusion, we have developed the first NN poten- 
tial at fourth order of xPT [40]. This potential is as 
quantitative as some so-called high-precision phenomeno- 
logical potentials. Due to its basis in xPT, the many- 
body forces associated with this two-body force are well- 
defined. Thus, we have a promising starting point for ex- 
act few-body calculations and microscopic nuclear struc- 
ture theory. 
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